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Abstract 

The measurability by means of continuous measurements, of an observable 
A(to), at an instant, and of a time averaged observable, A = 1/T J A{t')dt', 
is examined for linear and in particular for non-linear quantum mechanical 
systems. We argue that only when the exact (non-perturbative) solution is 
known, an exact measurement may be possible. A perturbative approach is 
shown to fail in the non-linear case for measurements with accuracy AA < 
AAmin{T), giving rise to a restriction on the accuracy. Thus, in order to 
prepare an initial pure state of a non-linear system, by means of a continuous 
measurement, the exact non-perturbative solution must be known. 
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I. INTRODUCTION 



The measurement of an observable A(to) is ideally described by an impulsive coupling 
between the system and the measuring device at t = to llll- In a realistic measurement 
however, due to the finite duration of the coupling between the systems, the back-reaction 
of the measuring device modifies the free evolution of the system. Thus, one generally 
observes a time averaged observable which depends of both the system and the measuring 
device. This raises the following questions: 

• Can the exact value of A (to) at the instant t = to, be measured by means of a contin- 
uous measurement? 

• Can the value of a time averaged observable, A = Jq A{t')dt', be precisely mea- 
sured? 

In this article we shall address these questions in the context of linear and in particular of 
non-linear quantum mechanical systems. Other aspects of continuous measurements have 
been discussed for example in Refs. 

The answer to the first question above has fundamental significance for quantum me- 
chanics. According to quantum mechanics, in order to prepare a certain initial pure state 
iplt = to), one has to determine the values of a complete set of observables at the instant 
t = to. But what if only continuous measurements are possible? 

It turns out that for linear as well as non-linear systems, as long as the exact solution to 
the Heisenberg equations of motion of the system is known, such a preparation is at least 
formally possible. To see this, consider a general exact solution for the observable A(t): 

A(t) =e(t-to,A^(to)), (1) 

which is given in terms of initial values of the observables A^ at time t = to. By sending 
t — s> — t and replacing the roles of t and to we can represent the constant of motion, A(t = to), 
in terms of time dependent operators Am(t): 

2 



A{t = to) = C{to-t,A„M- 



(2) 



We can now use the standard continuous measurement described by the interaction 



where Q is conjugate to the "output observable" P, and g(t) = ioi < t < to + T (or 
more generally, g has a finite support in time, and / gdt is finite). Since [A(to), Hj] = 0, the 
measurement does not cause an error in A{to), the exact value can be measured. However, 
this approach requires the knowledge of the exact non-perturbative solution. But what if 
the latter is unknown? We shall further discuss this issue in Sections 5-7. 

The second question raised above becomes particularly acute when asked in the context 
of quantum field theory. The singular nature of a field at a point requires to consider space- 
time weighted (smeared) fields as the elementary observables of quantum field theory 0. 
Therefore, we are faced with the problem of measuring a space-time averaged observable. 
The measurability of such objects has been investigated long ago by Bohr and Rosenfeld ^ 
for the case of a free field, but has not been considered for non-linear field theories 

As a toy model we shall consider in this article the problem in the simpler case of non- 
relativistic quantum theory. For the case of a linear theory, the back-reaction on the system 
during a measurement of A = T"^ / A{t')dt' depends only on variables of the measuring 
device. A compensation term can therefore be devised so as to exactly compensate for the 
back-reaction, and measure the latter averaged observable to any desired accuracy. But if 
the system has non-linear equations of motion, the problem of observing a time averaged 
quantity becomes much more difficult. The error due to the non-linear back- reaction becomes 
a function of both system and measuring device variables. As already described above, if a 
non-perturbative solution is known, then formally a precise measurement is possible. Since 
in general however this is not the case, we shall consider in some details the problem of 
measuring a time averaged observable by means of a perturbative approach. We shall 
show that the validity of a perturbative method (in the non-linear coupling constant) is 



Hi = -g{t)Q^{to-t,A^{t)) 



(3) 
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limited. For measurement with of accuracy better than a certain minimal uncertainty, the 
perturbation scheme breaks-down. 

In the next section we shall present the toy model which will be used in this article 
to study linear and non- linear harmonic oscillator with a non-linear potential. In 

Section 3. we shall consider three different approaches to measure averaged observables in 
the linear case. We then discuss these methods when the system is non-linear in Section 
4. A perturbative approach to the non-linear case is developed in Section 5, and shown to 
break-down for precise measurements in Section 6. We conclude by a discussion of the main 
results. 

II. NON-LINEAR HARMONIC OSCILLATOR AS A TOY PROBLEM 

We shall consider as a toy model for a non-linear system a harmonic oscillator with a 
potential given by 



Our first aim will be to describe a measurement of an averaged observable such as 



In the limit of T ^ the prescription is know: the appropriate interaction term is in this 
case 



where Q is canonically conjugate to the "output variable" P, and g(t) = go5{t). It is assumed 
that the effective mass of the device is infinitely large, thus the kinetic part of the measuring 
device Hamiltonian vanishes and Q is a constant of motion. The interaction Hamiltonian 
(|) yields: 




(4) 



where n > 2. 




(5) 



Hi = -^(t)Qx, 



(6) 
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5P = P(+e)-P(-e)=yox(0). (7) 

As a starting point let us modify this interaction by replacing the Dirac-delta function 
by a smooth function g{t) with a compact support only the time interval < t < T. 
The solution to the equations of motion 

x = p 

p = -Q2x + Ax"-^ + ^(i)Q 

P = g{t)^, (8) 

can be written in the integral form: 

x(t, g) = xo(t) + /* ^ sm[Q{t - t')]dt' (9) 

Jo Si 

p(^, g) = Po{t) + t F{t') cos[n{t - t')]dt', (10) 

Jo 

where 

F{t)^X^^-\t,g) + gm, (11) 

and Xo(t), Po(^) are free solutions {X = g = 0), which coincide with the non-linear solution 
x(i) ait — 0. The presence of the coupling parameter g in x(t, g) will be used in the following 
to denote an explicit dependence of the solution on variables of the measuring device. 

The solutions x(t, (/) and p{t,g) depend on the back-reaction via the term g{t)Q. Con- 
trary to classical mechanics, due to the uncertainty principle, we can not make Q as small 
as we wish and still obtain an accurate measurement, i.e. AP(0) — > 0. Consequently, we 
finally observe 

P(T)-P(0)= r g{t)^{t,g)dt, (12) 
Jo 

which is not the undisturbed (g=0) value of x. Only in the limiting case T — does the 
error vanish. 

In next section we show how the undisturbed x can be observed in the linear case; A = 0. 
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III. THE LINEAR CASE 



We shall now proceed to examine the case of a measurement of an averaged observable 
and show how to eliminate the back reaction in this special case. 

For simplicity we shall choose the weight function as g{t) — go for t e (0, T) and zero 
otherwise. In this case, the solution of the equations of motions when A = is 

, , , , 9o(l — cos fit) ^ 
x(i) = xo(i) + ^^^^—^^ 

= Xo(i) + Q;(i)Q. (13) 

Thus, the P coordinate of the measuring device will be shifted by 



= P(T) - P(0) = go[At')dt' = g, j\o{t')dt' + (l " (14) 



The undisturbed average x = Xq = T" Jq yio{t')dt' is given by 



6P go f sinl^T^. 



goT Vf^ V VtT 

Since Xq depends hnearly on both P and Q, a precise measurement of P causes a larger 
uncertainty in the second term above. I.e. it increases the back reaction of the measuring 
device on the oscillator. Since the uncertainty in x is 

the minimal uncertainty is 



Ax_ = ^^i-^^-4°y'^r, (17) 

which vanishes only in the impulsive limit. We also note that in the limit of T — > cxo the 
disturbance does not average out but rather approaches the constant ^. 

The direct approach therefore fails to measure Xq precisely. There are however ways to 
correct or eliminate the error above, and in the following we present three different ways to 
achieve this goal. 
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The idea of Bohr and Rosenfeld was to correct the error by adding to the Hamiltonian 
(^) a new "compensating" term. The error in the shift of P in eq. ( p!^ appears as hnear 
in Q, very much hke the effect of a spring in the Q-coordinates. Since the coefficient which 
multiphes Q in ( p!^ is known, it is straightforward to compensate for this error by adding 
to the interaction Hamiltonian (^) a spring term: 

Hj = -g{t)q^ + hQ\ (18) 



where k can be chosen for example as 



g\t)T / sinQT 



It is straightforward to see that the new spring-term in the equations of motion precisely 
eliminates the back-reaction of the device on the system. In the limit T — 0, the compen- 
sation is of course not needed since the spring constant vanishes, and we obtain back the 
ordinary impulsive measurement. 

Another approach due to Unruh [^, does not require modification of the interaction (P) 
but requires instead an additional measuring device. Inspecting eq. (^) we notice that xq 
is given on the right hand side in terms of a linear combination of P and Q. Therefore, after 
preparing the measuring device in an initial state with a well defined P one can measure by 
means of another measuring device the linear combination 

^P-^Q(T), (20) 

90-1- go 

where k is given in eq. (0). As before this method requires that the back-reaction effect on 
the system depends only on variables of the measuring device. As we shall see in the next 
section, both methods apparently fail when non-linearities introduce back-reaction terms 
which depend also on the system itself. 

A third different method to measure x could be the following. By integrating x()f:) we 
can express x as 

^ ^^(xo cos nt' + ^po sin nt')dt' = ^ sin fiT + ^(1 - cos ilT) . (21) 



Note that xq and po are constants of motion. By inverting the solutions for the equations 
of motion for x = x(xo, po) and p = p(xo, po) we obtain 



rx(t) = -smQT/2 

U Li 



x(t) COS n{T/2 -t) + ^p(t) sin n{T/2 - t) 



(22) 



Thus, we have expressed x in terms of a (time dependent) constant of motion. We can now 
set up a continuous measurement 



Hi = -^(t)Qx(t). 



(23) 



Clearly the observable x(t) remains a constant of motion and is not disturbed by the inter- 
action. We have thus managed to measure x by a continuous measurement of a constant of 
motion whose value is identical to x for any t. 



IV. NON-LINEAR CASE 

The presence of non-linear terms makes the back-reaction dependent on both system 
and measuring device variables. This makes the problem of constructing the analogous 
compensation more complicated. 

Nevertheless, if the exact solution is known, there is at least a formal way of eliminating 
completely the back-reaction. To this end we can use the third approach that was described 
in the last section. Let use denote the exact solutions of the equation of motion by 

x{t) = C(xo, Po, t), p{t) = ^'(xo, po, t), (24) 

Then, the average x may be expressed as 

Tx= re(xo,Po,tV^' = 0(xo,Po,T), (25) 
Jo 

where the latter observable, 0, is a constant of motion. We can now substitute 

xo = e(p(t), x(t), -t), Po = e'(p(t), x(t), -t), (26) 

into eq. (^51) and express in terms of the time dependent dynamical variables: 
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X = 0(e(p(t), x(t), -t), e'(p(t), x(t), -t), t) = 0(x(t), p(t), t, T). 



(27) 



Since 0(x(i), p(t), t, T) is a constant of motion we can measure it by means of the ordinary 
interaction (|^). 

The other two methods of Section 3. are much harder to use. Assume again the exact 
solution is given, and let us express it as 



Here, AC,(t, A, g) = C,(t, A, g) — C,(t, X, g = 0) = x{t, g) — x{t, g = 0) is the non-linear error due 
to the measurement. Thus, using the naive coupling we obtain: 



To proceed, we would like to construct a compensation to the last term.Q Thus, must be 
expressed in terms of the dynamical variables x(t, g), p(t, g) and Q. However, depends on 
x(t, (7 = 0) and we first need to find the non-linear relation between the disturbed non-linear 
x(t, g) and the undisturbed x(t, (7 = 0) solutions: 



etc. The latter equation corresponds to a "canonical" transformation between two non-linear 
solutions of two different systems. It is not clear if such a well defined relation indeed exist. 
Assuming however that transformation above does exist, we may now attempt, as before, 
to add to the Hamiltonian the compensation term /[A;Q + A,^(A, (7)](ig. However, the new 
equations of motions will give rise to another non-linear error A,^' ^ A^. Therefore, a self 
consistent scheme for constructing a compensation must be found, which takes deals with 
the non-perturbative effects of the compensation on the system. 



At this point we note that in principle we could have arranged that during the measurement 
a compensation term —g{t)Xx^ /n eliminates the non-linearity. We could then use the previous 
method for observing the time averaged observable corresponding to a a linear dynamics. 



x(t,^) = i{t,\,g 



0) + Ae(t,A,(7). 



(28) 




(29) 



x(t,^) = /[x(t,0),p(t,0),Q], 



(30) 
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The difficulties in finding an exact non-perturbative compensation scheme, leads us to 
examine a more modest approach. It can be hoped that at least for small averaging time T, 
the effect of the non-linearities is small. Therefore, in the next section, we shall examine a 
perturbative approach. 

V. A PERTURBATIVE APPROACH 

For a small non-linear constant. A, we shall attempt a perturbative approach. As a 
starting point we shall assume that the measurement can be expressed by: 

Hi^-g{t)QK + Hc, (31) 

i.e., as a sum of the naive measurement and a compensation term. Clearly in the limit of 
T — > 0, if/ should reduce to an ordinary impulsive measurement, i.e., Hc{T = 0) = 0. Our 
aim will be to construct a self-consistent procedure to compute the compensation term He 
to any order in A: 

He = AH^°^ + ^ X''AH^''\ (32) 

k 

The terms AH^''^ are the analogous compensations to the back-reaction up to the k'th order 
in A. 

To this end, let us expand the solution for x(t) as 

x(t, g) = x(°) + Ax(^) {t) + A'x(2) (t) + • • • (33) 

and at each order separate the undisturbed solution, denoted as x.'q \ from the non-linear 
error, denoted as Ax^''^: 

^^'\9)^^^\t)+A^(''\t,g), (34) 

i.e. for gr = we have Ax.^''\t, g — 0) — and x.^'^\g = 0) = ^'o\t) is the solution without 
the coupling to a measuring device. 
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Likewise, the shift of the output register , SP = P(T) — P(0), will be expanded in powers 
of A, and at each order we shall evaluate the terms of the decomposition: 

(^P^'^) = 5P('^ + AP('=). (35) 

Here, AP^'^) is the error due to the back reaction to the k'th order. The knowledge of AP^^ 
will allow us to construct the appropriate compensation /S.H^^\ up to the same order k. 

To proceed we shall now simplify the problem, and let n = 3 in eq. (^. It can be shown 
however that the following will be valid for any general non-linear potential. 

To evaluate AP^'^) we first use eqs. to obtain an integral solution to any order: 

xW(t,^7) = xo(t) + Q f g{t')D{t - t')dt\ (36) 

Jo 

x«(t,(7)= f\'^'>\t')^'^°\t')D{t~t')dt', (37) 
Jo 

^^^\t,g) = 2 f\^^\t')^^^\t')D{t-t')dt', (38) 
Jo 

etc., where xq is the free solution {g = \ = 0), and 

^, , sinQi 

m = (39) 

The result for the zeroth order was already obtained above, where we found 

AHo = hci^ (40) 

where k is the c- number given by eq. (0). Here, and in the following we shall ignore any 
potential problems with ordering of non-commuting operators. 
To the first order in A we obtain 

5p(i) =^0 Trft Ax(°)]2dt'. (41) 
Jo Jo 



Thus, 

2a(t')Qxo(t') + «'(i')Q' 
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AP(i) =go f dtf 
Jo Jo 



dt' 



^0 



2a(t')Qxo(t') + a\t')Q^ 



dt' r e{t-t')dt 
Jo 



{{T - t')e{T - 1') - {-t')e{-t') 

= go - t') 

Jo 

Jo 



2a(t')Qxo(t') + a'(OQ' 



dt' 



2a(t')Q(x^"^(t') - Qa(t')) + a'(OQ' 



dt'. 



dt' 



2a{t')Q^^^>{t')dt' - a\t){T)Cl^ 



(42) 



where a(t) is defined in eq. (|l^). As could have been anticipated, the new feature of the 
error terms up to the first order in A, is their depends on the variable x of the system. The 
last term does not depend on system variables and therefore may be trivially compensated as 
a spring-like compensation. (Alternatively, as discussed in Section 3., we may use a second 
measuring device to measure the combination — kQ — A;'(T)Q^, where k' can be found 
by integrating the last term in eq. (^21). ) 

Equation (^) suggests adding as a compensation the term 

1 



Ai/« =/(t)(r-t) 



a(t)QMt) - 3«'(t)Q' 



(43) 



With the compensation AH^^'^ + \AH^^^ added to the Hamiltonian, the eqs. of motion 
become 



P = (?(t)x-A;(t)Q + A/(t)(r-t) 



2a(t)Qx(t) -a2(t)Q^ 



'd^ 



+ n^x = Ax^ + g{t)Q + W{t){T - t)a{t)Cf. 



(44) 



(45) 



The new term on the right hand side of (|45D modifies the solution for x{t) to 

x(t) = xo(t) +a(t)Q + A/5(t)Q2 + A ^^{t')dt'. 

Jo 



(46) 



The first order is modified by the /3(t)Q^ which can trivially be compensated by adding also 
the term |A(7(t)/5(t)Q^. This completes the compensation to the first order. It is important 
to note that the a first order compensation can not modify lower order compensations but 
only higher orders. This allows us to go to higher orders in A without effecting the corrections 
found at lower orders. 
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Let us see how one can proceed to higher orders. To the second order we obtain: 



dt 



where 



g{t)QA^^'^ + 2g\t){T - t)a{t)ClA^^^>{t), (47) 



Ax(^) (t) = + f A[x(°) it')?dt', (48) 

Jo 



and 



Ax(2) (t) = 2 f A[x(°) (t')x(^) it')]dt' (49) 

where the notation A[---], on the right hand side in the eqs. above, indicates that only 
Q-dependent terms of [• • •] are included. 

The computation can be carried out, as before. The only technical subtlety is that now 
AP*^^-* contain also triple integrals over time like: / dt J XQ[t')dt' J pf{t")dt" . To proceed 
one has to reduce this integral to a single integration over time, or eliminate the integrals 
completely. Otherwise, the compensations would have to be non-local in time. In our 
case of a first quantized theory, we can always use the property of linear operators (but 
not field operators): po{t) (or xo(t)) can always be expressed as a linear combination of 
Xo(i = 0) and po{t — 0). Therefore, we can integrate over free operators and reduce: 
/o Po{i")dt" = ci(t')xo(t') + C2(t')po(t')- The integral above can be hence reduced to a single 
integral of the form J{c[{t)xl + C2(i)xoPo)c?i. Finally we find that the compensation to the 
second order has the following form: 

where 7i(T) are generally time dependent c-numbers. In principle it seems that this proce- 
dure can be carried out to any order in A. 

VI. BREAK-DOWN OF THE PERTURBATIVE APPROACH 



Although by using a perturbative approach we can formally evaluate self-consistently 
compensation terms to any order in the nonlinear coupling constant, this approach still does 
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not allow precise measurements. Unlike classical measurements, in the limit of a precise 
measurement, the variable Q, which is conjugate to the "output register" P, becomes com- 
pletely uncertain. Since the expansion for the compensation He generally depends on Q, 
for a given T, the validity of the perturbative expansion must break-down beyond a certain 
precision Ax^m- 

To show this consider the first order compensation term. From eq. (^3|) we have AH^^^ ^ 



Ta(T)Q^x. A necessary (but not sufficient) condition for a well behaved expansion is 
therefore: 

ATa(T)Q2x ~ AT^^Q^x « 1. (51) 

Substituting a rough approximation ~ (^Q)^ and x ~ Ax and using AP ~ Ax/c/o^^ we 
obtain the condition 

Ax > A^min ^ XglT^ (52) 
Only by letting T — > or A ^ we regain Ax^m — > 0. 



VII. DISCUSSION 

In this article we have examined the measurability of observables by means of a continu- 
ous measurement for linear and non-linear theories. We have seen that as long as the exact 
non-perturbative solution is known it is possible, at least in principle, to device an exact 
measurement. The formal answer to the two questions asked in the introduction is therefore 
positive. 

Nevertheless, the surprise is that without the knowledge of the exact non-perturbative 
solution, a perturbative approach can not be used to obtain arbitrary accurate measure- 
ments. For any non- vanishing value of the non-linear coupling constant, and a given in- 
teraction time T, the perturbation breaks-down for any measurement with accuracy better 
than A A < AAmj„(T, A). The origin of this failure is that the perturbation to any or- 
der, necessarily involves terms which depend on Q, the conjugate to the output variable 
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p. Therefore for sufficiently precise measurements, AP, can be arbitrarily small, and AQ 
becomes arbitrarily large. At a certain accuracy the perturbation is bound to break-down. 

This conclusions clearly applies to the question of measurability of an observable at an 
instant by means of a continuous measurement. If the exact solution is not known, we shall 
not be able to compensate perturbatively for the error. Thus, under the restriction of finite 
time measurements, a preparation of a certain pure state if) oi & non-linear system, requires 
the knowledge of the exact non-perturbative solution of the equations of motion. 
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